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Abstract 


‘¥e~~e*4umd earlier plasraajklnetic model of an inverted-V auroral arc 
structure ^ include^i-n-«~phenQmenologle,al'^ji^ the effects of electrostatic 
turbulence ^ In the absence of turbulence, a parallel potential 

drop is supported by magnetic mirror forces and charge quasi^eutraliUy, with 
energetic auroral ions penetrating to low altitudes; relative to the elec- 
trons, the ions’ pitch^^angle distribution is skewed toward smaller pitch 
angles/ 'fhe electrons energized by the potential drop form a current which 
excites electrostatic turbulence; — we--^onsidert>4;hQ-*spec4'f-i^=ca3er''H3f~tho= jfon*- 
Gy-cl-c>OroTr"m^aV .Jaifis=taWifl«nee^lte^S~"tTle^ thus<=tending 

redyjse^ the-„d±£fere.atial pit.di , angle anlaotrorpy ’beliween electrons and ionsf^s; 
whichL.ln tura reduces the potential dro^ - .an-eMeet of., opp osite sign to_ .ttgcsy**" 
a-ssCTcl'ared with anomalous reslstlvityv^ In equilibrium the plasma is marginal- 
ly stablq^ wii±„giOLWiJbL.rntes, a diffusion cons tanta-sorae—two-oi-ders-of magn^P^^- 
(?Ude— belW nai-ve ..a st ima tsaaa-.. The conventional anomalous resistivity contribu- 
tion to the potential drop is very small* ■‘bec«tt^^e _anomalous^esistivity pro- 
cesses are far too dissipative to be powered by auroral particles ;-4A4s— 
->-g^?owth“Trartes~'^TTrffu:s±on--Gonstaiita,-nre ..sO'-smaS^ Under certain circum- 
stances equilibrium may be impossible and relaxation oscillations set in 
t-irae' scale -foT^ffch“pulsatdni:g^-aurar as is the ion t ransit tim e„Ouf--A=U^--BeeeSda^l> 


Abstract 


We extend an earlier plasma-kinetic model of an inverted-V auroral arc 
structure to include, in a phenomenological way, the effects of electrostatic 
turbulence with 1* In the absence of turbulence, a parallel potential 
drop is supported by magnetic mirror forces and charge quasi-neutrality, with 
energetic auroral ions penetrating to low altitudes; relative to the elec- 
trons, the ions' pitch -angle distribution is skewed toward smaller pitch 
angles. The electirons energized by the potential drop form a current which 
excites electrostatic turbulence; we consider the specific case of the ion 
cyclotron mode. The turbulence heats the ions in Tj^ only, thus tending to 
reduce the differential pitch angle anisotropy between electrons and ions, 
which in turn reduces the potential drop - an effect of opposite sign to that 
associated with anomalous resistivity. In equilibrium the plasma is marginal- 
ly stable, with growth rates and diffusion constants some two orders of magni- 
tude below naive estimates. The conventional anomalous resistivity contribu- 
tion to the potential drop is very small, because anomalous-resistivity pro- 
cesses are far too dissipative to be powered by auroral particles ; this is why 
growth rates and diffusion constants are so small. Under certain circum- 
stances equilibrium may be impossible and relaxation oscillations set in; the 
time scale for such pulsating auroras is the ion transit time of 6-10 seconds. 
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I. INTRODUCTION 

Recently many authors have been concerned with plasma-kinetic mechanisms 
for explaining auroras associated with inverted-V's (auroral arcs, for brevi- 
ty) [Swift, 1975; Kan, 1975; Hudson and Mozer, 1978; Chiu and Schulz, 1978; 
Fridman and Lemaire, 1980; Lyons, 1980, Chiu and Cornwall, 1980]. The funda- 
mental objective is to set up a parallel potential drop of 1-10 keV along an 
auroral field line which will accelerate electrons and precipitate them into 
the ionosphere. To this end, three mechanisms are most widely cited: 1) 
charge-separation effects (Poisson's equation); 2) anomalous resistivity; 3) 
magnetic mirror forces coupled with differential electron-ion pitch-angle 
anisotropy. Because of the large (> 10^) transverse dielectric constant of 
the magnetized auroral plasma, transverse derivatives in Poisson's equation 
completely dominate derivatives along the field line [Chiu and Cornwall, 
1980]. We will therefore not consider unmagnetized double layers (potential 
drop ~ keV in a distance ~ 0.1 km along a field line) any further, although 
several authors have invoked them for auroral arcs [Block, 1975; Shawhan et 
al., 1978]. The remaining effect of Poisson's equation is to couple neighbor- 
ing field lines on a length scale R^ (Rl is the ion Larmor radius) [Swift, 
1975; Chiu and Cornwall, 1980]. This sort of coupling will not be vital for 
our present purposes, so we will demand instead that the plasma be quasi- 
neutral: Nj^ =» Ng everywhere. 

Next consider anomalous resistivity. While it may be very important in 
some laboratory circumstances, a simple energy-balance argument severely 
limits its role in the aurora, at least in the scenario we adopt. In that 
scenario, all auroral energy comes from the kinetic energy of plasmasheet 
particles, and is likely to be dominated by the ions. Suppose that turbulence 
somehow generates an anomalous collision frequency v. for electrons (of mass 
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in) I so Chat the curreat-fielcl relation is 

2 

J| " 0|| E|| } 0|| " Ne /m (1) 

The electron current is of order 4, « where is a typical electron 

velocity, so that the power dissipated on the auroral field line, integrated 
over the region of length A whore the turbulence is present is, in order of 
magnitude, 


/ ds 4j/0|| « A Nm (2) 

This must be less than the input power, which comes from fields and particles 
in the plasmasheet. There are several estimates one can make for this, but 
none are plnusibly larger chan the ion energy flux. (Quasi-neutrality of the 
plasmasheet limits the electron energy flux to be of the order of the ion 
energy flux.) For example, one might begin with the power per unit volume 
* ^1 dissipated in MUD flow. We estimate 4^^ cVp/B , Ej^ ~ 1 mV/m, and 
multiply the product by a length which we take to be the same as that setting 
the scale of the press are gradient. The result is substantially less than the 
ion power NMV^ (M is the ion mass; V is a typical ion velocity), and it is 
this latter quantity which we use as the input power; it is of the order of a 
few ergs/cm sec. In terms of the average plasmasheet electron energy and 
plasmasheet ion energy W, we have >• mW/(MWg), so that requiring (2) to 

be less than NMV implies 

Va < (W/W^)(V/A) (3) 
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For typical auroral panunetora this amounta to < 10 , considerably 
less than the range " lO"^ - 10“^ found in simulations or laboratory 
experiments, which refer to quite different physical circumstances not found 
in auroras. 

The point, however, is not the absolute size of v^, but rather what sort 
of potential drop it can lead to. Anomalous resistivity is so dissipative 
that if turbulence and heating roach the maximum level allowed by (2), only a 
small potential drop is generated (a somewhat similar point was made by 
Falthamraar [1977]). Let us use (2) and (3) in (1), estimating the total 
potential drop between the ionosphere and the equator as A(j> “ -Z K|| . One 
finds 


OM < w 


1/2 


(4) 


In the plasmasheet, W/Wg " 2-4 is typical [e.g., Frank et al., 1981] so the 
potential drop, scaled in units of ion energy, is perhaps 0.05 - 0.1. Note 
that this is independent of the length J of the turbulent region. In contrast 
to the bound (4) on eAij) maintainable by anomalous resistivity, the theory 
outlined below, as well as experimental observatlone, suggest that eA<j) " W. 


Thus the anomalous-resistivity potential drop is limited severely by the 

available auroral power, which is why canonical efitimates " (10““^ - 

10 ) 0 ) are so far off the mark: they refer co circumstances where the 

pc. 

potential drop is maintained by an external battery, which can supply the 
power which is necessarily dissipated by anomalous resistivity. 


This is an example of a theme which will occur several times in the 
course of this work; microscopic wave-plasma processes are capable of trans- 
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fenring energy between particle epccios, or turning this energy Into wave 
turbulence, at rates which are enormous compared to the Input power avail- 
able. The plasma processes are analogous to an enormous pipeline through 
which only a trickle of water flows; it is completely wrong to estimate the 
water flow from the diameter of the pipe. Tlie overall balance of plasma 
processes is in large part governed by macroscopic, not microscopic, para- 
meters. (This is true at least for plasma effects in which > 1, but not 
necessarily if kj^ < 1.) 

Turn now to the third kinetic process of magnetic mirror forces. It has 
been considered in great detail by Chiu and Schulz [197SJ (using 'luasi- 
neutrality) and Chiu and Cornwall [1980] (who added Poisson's aquation and 
ionospheric charge and current conservation, effects whose greatest importance 
is to couple neighboring field lines). These authors (as well as Lyons [1980] 
and Fridman and Leraaire [1980]) pointed out that there is, for magnetic fields 
with large mirror ratios, a linear relation between parallel current and 
potential drop 


- ( 5 ) 

(we drop a term referring to diffuse particle precipitation) with Q " Ne^/(m 
Vg). Again by taking A(|) « “f-Ej , one finds by comparing (5) with (1) that the 
effective "collision" frequency is simply the inertial frequency 

Vj » V^/Z (6) 

which is larger by roughly the factor (M/m)^’/^ than in (3), and of course 

2 

is consistent with eA(|i " MV , as observed. It is important to note that this 
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much lacgcr effective col. Ision frequency 1 b not associatea with diooipation 
of any sort, and equations such as (2) with Vj. in place of are meaningless 
for the magnetic mirror force process. 

Somewhere between the extremes of totally dissipation-free auroral arc 
models (e.g., Chiu and Schulz [1978]) and moximally-dissipativc models of 
anomalous resistivity lies a middle ground in which there is turbulence and 
dissipation, but in an amount limited by the auroral input power. We have 
already argued that this limitation implies little or no turbulence contribu- 
tion to the parallel potential drop, so we must look elsewhere for the effect 
of turbulence. In fact, it turns out that the effect is precisely opposite to 
that of anomalous resistivity: electrostatic ion cyclotron turbulence tends 
to reduce the parallel potential drop. The argument is that such turbulence 
(or similar electrostatic modes with E| « is very efficient at heating 
ions, thereby increasing Tj^/Tj for the ions. (At the same time electrons are 
slov;ed down and heated mostly in T^). As a result, the differential pitch- 
angle anisotropy betv/een ions and electrons is reduced. In the Chiu-Schulz 
kinetic model, precisely this differential anisotropy drives the parallel 
potential drop, which - as the detailed calculations presented below show - is 
reduced as T,/T„ increases. 

1 n 

There are now two effects tending to stabilize the turbulent wave ampli- 
tude: 1) the well known (e.g., Palmadesso et al. [1974]) Increase in the 
electron critical drift velocity with increasing T^^; 2) a decrease in the 
electron drift velocity Vy because tVie parallel potential drop has 
decreased. As Palmadesso et al. [1974] show, the first effect alone acts to 
bring Vq - to a very small value (compared to other characteristic veloci- 
ties), and the ion cyclotron mode is essentially marginally stable (the growth 
rate Y is linear in Vq - Vq). With only the first effect present, we will 
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show below an exponential decay to marginal stability (Y“0) » Adding the 
second effect may have more dramatic effects s with Vy decreasing at the same 
time Vq increases it is possible for y to go negative, with consequent shut- 
down and decay of turbulence* Eventually fresh ions (with ^ 

injected and restore the initially unstable situation. These ions enter the 
system on the ion transit time scale Ji/V » 10 sec; since this is greater than 
the time scale for shutting off turbulence, the arc can pulsate at this peri- 
od. In the present paper we will not take up this mechanism for pulsating 
arcs in any detail} however, we note that the observed characteristic period 
of pulsations (both patches and arcs) is very nearly the ion transit time 
[Royrvik and Davis, 1977]. 

Whether pulsation or decay to marginal stability takes place depends on 
Initial conditions} eventually some sort of equilibrium will be reached. We 
can make some simple estimates of turbulence parameters for the electrostatic 
ion cyclotron mode, when a marginally stable equilibrium has set in. First it 
is necessary to identify the wave saturation mechanism. Formation of a quasl- 
linear plateau on the electron distribution is not favored, because old elec- 
trons are rapidly transported out of the region of instability while now ones 
are being carried in, and convective losses are slow because of the very low 
wave group velocity. It is thus plausible that ion resonance broadening [Dum 
and Dupree, 1970] and the rather smaller effect of geometric resonance broa- 
dening (e.g., SchulH [1972]) are the dominant saturation mechanisms. Dum and 
Dupree show that, in this case, the RMS wave amplitude is largely independent 
of the ion velocity diffusion coefficient} their formulas applied to the 
auroral arc problem suggest vmve amplitudes of the order of tens of mV/m, 
which is roughly what is observed [Temerin et al., 1981]. The diffusion 
coefficient for ion diffusion in Vj is related not to the wave intensities, 
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but to the linear gcowth rote 


where fl is 
- fi/v, 


n * 

1) a iL^ Y 


the ion cyclotron frequency# The most unstable modes 
and Yt is of the order of 

14 


Y, . Tja (V, - v,)/v„ 


(7) 

have 


( 8 ) 


where ^ 0# 2, Vp is the electron current velocity along the field line, end 
Vq, the critical drift velocity, is nominally of order 10 V « 0.2 [Kindel 
and Kennel, 1971). Thus a nominal value for Yj^ is of order 0.1 rj; with nomi- 

nal auroral values of 0 “ 100 sec ^ and V >=* 10° cm/sec, wo estimate D » 10 

3 -1 

cm*' sec 

In fact this estimate for D is too large by orders of magnitude, for 
essentially the same reason as the nominal estimates of anomalous resistivity 
were wrong: such a D yields too much dissipation to be supplied by the a*,sro- 
ral input. As an ion travels along the field line from equator to ionosphere, 
its perpendicular energy increases by wave turbulence, by an amount roughly 
equal to 2M j ds D/V|j . Since the ultimate source of energy for ion heating 
is Che ions themselves (that is, heating in is a less-than-perfectly 
efficient way of transferring energy from parallel ion motion to perpendicular 
ion motion, via a potential drop which induces turbulence) this quantity must 
certainly be less than roughly the ion kinetic energy which gives the limit 

D < ^ » 10^^ cm-^ sec ^ (9) 
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(oDCimatins « 4 R^). When thlo upper Umlc aa D ia uoed in (7) to oatimate 
Yj^, it ia aeon that Yj^ oeaXoa ao V/X, « 0.1 it. Thio meano, of couroe, that 
turbulence io generated at very near marginal stability. 

In this paper we make quantitative the considerations given above, using 
a ochematic but adequate representation of electrostatic turbulence with 
kj^ » k|. The starting point is the Boltcmann equation for particles in an 
inhomogeneous magnetic field, with on applied 00 potential drop, and In the 
prooenco of perpendicular velocity'-opace diffusion. The diffusion coefficient 
D is Independent of velocity, and is considered as an appropriate average over 
the turbulent spectrum. Addition of turbulent diffusion makes the Uoltaraann 
equation analytically insoluble, but we show that for an initially bl- 
Maxwellian velocity distribution function f, a good approximation is that f 
continues to bo flaxwellian, but with a temperature T., as well as an of foe- 
tive electrostatic potential, dependent on field-^^llne incegrals of I) weighted 
with various functions. (T^ does not change because of D.) That is,, the 
effects of diffusion c«n bo well modeled by leaving out the explicit diffusion 
operator and renormalixlng T, and the parallel potential drop 

This is a great simplification, because it allows us to treat our problem 
as a simple modification of the Chiu-Schulz problem (to the extent, at least, 
that we ignore coupling of neighboring field lines and the ionosphere, which 
is permissible here). We have run the Chlu-Schulz numerical programs for 
various values of D (that is, various renormalized Tj^ and A(j>) and find that, 
as anticipated, A(ji is forced to decrease as D increases. At a critical value 
of D the solutions either break down or Aif is essentially zero; the critical 
value of D found in the numerical work is somewhat smaller than given in 
(9), This is, of course, expected since (9) is an upper limit. 
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Section II contains a derivation of the approximate solution of the 
Uoltsmann equation. Section III dioeuoooo the physics of oloetroototie ion 
cyclotron fu'buleneci and Section IV contains numerical results. Section V lo 
a brief recapitulation and statement of work for the future. 
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ll, BOLTZMANN liUUATION WITH YURBUUiUyE IN AH INHOtlUGlvUEOUB MAUNIvYlC FIELD 


Tho Boltamann equation for lono ouf ferine voXocity“npnee diffunlon in the 
component of velocity pcrpondxcular to the magnetic field io, after cyclotron 
averaging, 


0 m 1^ V i— + V 
^ l3t i an 


Fa 


I 3V. 


+ V, 


J_ 

9V. 


L 

Vi 


9V, 


(DV 


1 3V 


’)} f 


( 10 ) 


Where 





M + ii 

9fj M 


V- V 


V, « 


3B 

2B 3 s' 


( i ia ) 
(lib) 


We take D to be independent of velocity and time, and concern ourselves only 
with static (3/St ^ U) cases. For D « 0, the well known solution to (10) is 


‘ “ *0 "j,0> 


®tio 


M M V,^' B 

(1 - + e (<}. - <|.y) 


^10 



its 


2 


B 

B 


0 


( 12 ) 


(13) 


(14) 


The subscript 0 labels a particular point on the field line, where s « 0, 
B « Bq, and the electrostatic potential “ ifiy (see Fig. 1). In the case of 
energetic plasmasheet ions, the point 0 will be the equator of an auroral 
field line. 
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If U 0 it is impossible to solve (10) exactly, except in a formal 
sense. The formal solution is nonetheless instructive for further develop- 
ments, so we give it hero. 

First introduce the fundamental quantity 

U,(s, s', V,, Vj) « {vj + v'f (1 -|^) +|S. (4. - (15) 


where B' “ B(s'), B “ B(s) and analogously for ^ and V, and V|| , Vj^ ore the 
velocity components at s. We restrict ourselves to regions of space and 
velocity space such that Ujj is strictly positive, that is, no ion mirrors in 
the region between s and s'. The general case without this restriction re- 
quires very cumbersome notation. In our coordinate system (Fig. 1), > 0 
refers to downgoing ions, and in the earth's field we take s > s', B > B'. 

Clearly, Uj| is a constant of the motion for D « 0: 


L U 


In terms of Uj|, 


2 f V .L_ a. V q. V n ‘^OrUfs'asi^V 

'>11 8s i 8V, II 8V/ II ’ a''® I 

X a 

we define the fundamental diffusion variable Q: 

A W II 


(16) 


(17) 


and note that 


L Q = D (18) 


and in terms of Q, we define a Green's function for the diffusion operator: 
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( 19 ) 


V. VI 


“■• 3)^0 i 4 r ) ^ <''1 + 


G* (Vj^-vp, Q-0 

Ju 


( 20 ) 


Here 1 q is a Bessel function of imaginary argument. G obeys the Boltzmann 
equation (10) up to an error which is second order in s, the length of field 
lino over which the solution applies. That is, 


''i ’''j. 


(dv. li-) 

1 3V,'' 


•s 0 


[D (q"^ 1^-)^' Gj 


( 21 ) 


with the RHS of (21) failing to vanish because 3Q/9Vj^ 0. But from the 
definition (15) of , it is clear that, in order of magnitude, and for s << 
B/OB/36) 5 Iq, 


Q 


■1 3Q_ 
3V. 


U 


-1 


3U, 


n 3V, 


VJl, 


( 22 ) 


so the RHS of (21) is of order DG(s/VAq)^. A typical term on the LHS is 
0 (VG/Aq), and (as will be clarified later; see Eq, (9)) since D/V^ < V/A^, 
the error in the Green's function solution is (s/Aq) . Choose s to make 
this as small as desired; then, to the chosen accuracy, the solution to the 
Boltzmann equation which approaches a given boundary value fQ at s ■= 0 is 


f = / 


0 


n “''I ''P 


M V’^ B„ M V.'^ li. 

>“ £n 1-y- + ^ (1 - ir) + (t-to). -sJ 


M V 


BJ 


(23) 
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One may now simply iterate (23) going from s to 2s, using f at s as the input 
on the RHS, and then to 3s, ••• Ns. The exact solution is approached as s + 

0 and N + » with Ns fixed. 

This formal solution is of little use in general, but one special case 
will draw our attention; if £q in (23) is bi-Maxwellian, the output distribu- 
tion function is also bi-Maxwellian, except for the velocity dependence of 
Q. It will turn out that this velocity dependence is relatively unimportant, 
so we consider the possibility of an appruxiraately bi-Maxwellian solution by 
direct substitution in the Boltzmann equation. The solution ansatz necessari- 
ly involves an overall factor dependent on s, which can be interpreted as a 
sort of renormalization of the electrostatic potential. An alternative devel- 
opment of an approximate solution to (10) leading to similar solutions below 
can be obtained by the method of moments. Since this method may provide 
alternative insight into the solution scheme, we give a summary of it in the 
Appendix. 

The form chosen for f is 




(24) 


where 0|| may depend on s, but not on V. When this is substituted in (10) 
(with 8/9t = 0) the following equations are required to be satisfied, in 
addition to 0|| = constant. 


L0^+V,(6i-e,) 


1 9B , , „ 
B « 



0 


(25) 


Ln + 4 0j^D = O (26) 
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where h Is the differential operator defined in (i6). For U ° 0, the solu- 
tions are n " constant, 0„ constant, and 


Oj^ - 0j_(s) 5 (0j_Q - 0,) By/B + 0^ 


(27) 


Of course, this simply reproduces (12) - (14) for a bi-Maxwellian. 

For D 5^ 0, (25) and (26) are not strictly consistent, because of the non- 
trivial appearance of V^. To achieve consistency, we will replace. Vj| by an 
average value Uj which is independent of the local velocities V|| , Vj^, but 
might depend on s. In order that this approximation make sense, we must 
separately consider downgoing (V|j > 0) and upgoing ions; for the moment we 
treat only downgoing ions by taking f to vanish for < 0. Then U||(s) can 
be defined by averaging over this distribution function at various points 
along the field line, using the D «=> 0 solution for f as given by (12) - (14). 

In this approximation, Ou is still constant and (25) and (26) become 


3 

3s 


( 0 , - 


0,,) 


JL IM. 

B 3s 


+ 4 


2 - -1 


- 0 


(28) 


+ 4 0. D U„ 
3s 1 n 


(29) 


These equations appear to be non-linear, but actually they can be linear- 
ized. The first step is to multiply (29) by 0j^ and subtract it from (28); the 
result can be rearranged to 


3 

3s 


(B 


e-") 


-n ^ 
3s 


which allows us to express 0j^ in terms of Y = e 


( 30 ) 
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e,(s) 


y"^s) 


[0, / " ds' Y(s') H, + B 




(31) 


Now (29) can be written 


H 

9s 


£3 


4 U 0j_ U 


-1 


Y 


(32) 


and the combination of (31) and (32) gives a linear second-order differential 
equation for Y: 


as ^ 4 D 


3Y 

8s 


> “ ®i 



(33) 


This could readily be solved numerically, but the approximations made i.-' 

deriving the equation do not Justify anything more elaborate than saving terms 
of first order in D. (The reason is that we have assumed 0j^ and n to be 
independent of V. If this assumption is dropped, (25) and (26) could be 
solved without the approximation V|| U||. But then the equations themselves 

would not be quite correct; there would be terras of O(D^) coming from the 
action of the diffusion operator on 0j^ and n)» 

The lowest-order solution for 0j^ is 

° (4 0^5,-^ (s') 13^0^+®! <34) 

where, as before, B = B(s) and B' = B(s'). It is useful to phrase the lowest- 

order correction to n in terms of a phantom electrostatic potential <t)jj 

defined by 
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(35) 


‘<’0 


a (}) - 


Mn 

2o0, 


We find 


‘*’d " ° Ir / ^ ds' (D U^”^) (s') (30) 

II 0 

The RUS of (34) and (35) shows that tlie smallness parameter of our appro- 
ximation is roughly / ds' D0j^/U|| which is essentially the ratio of the 
perpendicular energy gain of an ion due to diffusion in a single pass along 
the field line (as noted in Sec. I ) to the initial ion energy. As we have 
seen, this cannot be larger than one, or the ions cannot maintain the necessa- 
ry differential pitch-angle anisotropy which drives the potential drop whicli 
drives the electron current which drives the turbulence. On the other hand, 
the smallness parameter is not necessarily much less than one, so our efforts 
here can be considered only seraiquantitanive. 

In summary, the primary effect of diffusion xb to leave a bi-Maxwelllan a 
bi-Maxwellian, but with modified 0j^ and with an effective electrostatic poten- 
tial <|)jj which is greater than the true potential drop. The Inverse tempera- 
ture decreases as diffusion acts, that is, Tj^ is increased as we expect. 
The fact that is greater than (fi might be interpreted as the presence of 
positive phantom charge which tends to keep the protons out of the ionosphere. 

The above discussion also holds with but trivial changes for upgolng 
ions, whose T^^ increases as they move toward the equator. Whether up- or 
downgoing, Tj^ heating raises the ion mirror point and diminishes the parallel 
potential drop; thus turbulence acts with opposite sign to anomalous resistiv- 
ity, and furthermore is self-limited by global considerations, not by the 
usual local saturation mechanisms. 
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To be complete we ought to look at the Boltzmann equation for elec- 
trons. As discussed in the next section, electrons are mostly heated in 
Tj| as well as having their field-aligned current velocity slowed down; that 
is, the flow motion is thernalized. This acts also to reduce the ion-electron 
differential pitch-angle anisotropy and to diminish the parallel potential 
drop. In view of some uncertainties in the quantitative description of this 
process (see Section III), wo do not attempt to treat electron heating and 
slowing down here, but it is certainly no more complicated than what we have 
done for the ions. We will be content to note that the effects of turbulence 
on electrons reinforce those acting on ions . 
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m. GLOBAL PLASMA PHYSICS OP THH ELECTROSTATIC ION CYCLOTRON MODE 


So far wo have only diacuaseii one-half of the problem; what effect dooa 
a given perpendicular- diffuaion conatant D have on the Ion diatrlbution? The 
other half la, of courae, how much turbulence and vduit value of I) do the 
diatrlbution functiona, magnetic field geometry, and boundary conditiona 
yield? 

It turns out that many of the plaama-phyaical numbers are largely deter- 
mined by Rlobal scaling laws, having little to do witli tlio local plasma envi- 
ronment. For example, the net growth rate y scales as the inertial frequen- 
cy V/A which involves A, the length of the field line, and D is of 

3 

order V /!L, The dominance of global over local effects will be an interesting 
challenge to experts in computer simulations of plasmas (recant work on ion 
cyclotron turbulence is reported by Okuda, Chang and Lee [1901]; Okuda and 
Ashour-Abdalla [19S1]). 

Let us briefly recapitulate eur scenario- Incoming auroral ions have 
average mirror points which would be closer to the ionosphere than those of 
the electrons, if it were not for the electric field which this enforced 
charge separation would produce. The parallel potential drop causes the 
electrons to have a net flow along the field line, which is strong enough to 
trigger the electrostatic ion cyclotron mode. The resulting turbulence (for 
which « 1) heats the ions in Tj^, which raises their mirror points and 
reduces the potential drop, thus reducing the electron drift velocity and. the 
growth rate. At the same time, ion heating in Tj^ acts directly to reduce tlie 
growth rate by increasing the threshold drift velocity (an effect invoked by 
Palmadesso et al. [1974]). One of two effects ensues; 1) the plasma reaches 
marginal stability in less than the ion inertia time i/V, with final values 
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of Y and U as estimated above; 2) the election drift velocity decreases and 
the critical drift velocity increases in such a way that becomes 

negative, the waves arc decayed or convected away, and the turbulence disap- 
pears. The initial unstable state will be restored after a time of order 
A/V, and the whole structure will undergo relaxation oscillations (observed 
as pulsations). 

In this section, wo discuss the plasma physics of the electrostatic ion 
cyclotron mode , including saturation mechanisms (ion resonance broadening, 
geometric resonance broadening); heating and slowing-down rates for ions and 
electrons; influence of this heating on growth rates; and dynamics of the 
approach eltlier to marginal stab**ity or to a relaxation-oscillation mode 
(pulsations). Let us begin with the dispersion relation for drifting bi- 
Max we 11 Ian distribution functions, thus generalizing Kindel and Kennel 
[1971 J. As in earlier sections, electron quantities are subscripted, and 
quantities without subscripts refer to ions, except in aquations (37) - (38). 

In a spatially homogeneous magnetic field the electrostatic dispersion 
relation is 


k « - f 


r.,(p) [u + Nfl(l- T.T ,“^)] - , 


rsa 


k. v„ 


k„ V„ 


(37) 


The sum is over all N from - to », and also over species. For each species, 
we define 


k 2 - 2 
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2 m 
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(38) 
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We take it that for ions there io no drift, and that (w + NO) » kj Vj , while 

for the drifting electrons w + NO « kj Vg . The fastest-growing modes have 

2 2 

n •< 1 (for ions), and kj /kj^ « 1. Then the dispersion relation is 

“2. a ( I (39) 

JL 


” T, X, U - - h (i - Tq)] 
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(40) 
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k V ^ 
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“ i)j 


(41) 


where p “ (u-0)/(kg Vg). The linear growth rate in a homogeneous field is 


thus 


Y C3 


X 0 Tg 



(42) 


In an inhomogeneous field, there is damping from geometric resonance broaden- 
ing, that is, ion cyclotron resonance is lost as the wave propagates along the 
field. A straightforward calculation shows that the damping decrement is 





1/2 

Hi 

8s-' 


(43) 


this is of order 1 sec"*^ on auroral field lines, and is not insignificant. 
The net growth rate is y ” ■" 

For this mode, the group velccity is small compared to the electron drift 
speed, so that fresh electrons are constantly being injected into the region 
of instability. This inhibits quasi-linear saturation by formation of a 
plateau on the electron distribution function, so the most important satura- 
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tion mochnnioroo arc ton reaonaneo broadening [Piun and Puprcci 1970] and ion 
heating [I’alraadeaso ct al., 197A). According to Dura and Dupree, the damping 
decreraont duo to ion rooonanco broadening io 


Y 


a 

iiib 




(44) 


and equilibrium is reached when Y “ Y^^y* hero D is the diffusion coefficient 
for ion diffusion in as used in Section II. As mentioned in Section I, 
ostimates such as Vy - u/k| « m/k|, or Y ~ 0.1 0, yield estimates of U which 
are much too large to bo furnished by auroral ions. 

In addition to ion diffusion (heating in Tj^), the electrons suffer a 
slowing-down of their drift velocity and heating in Tj|^. These processes are 
non-resonant (see the inequalities below (38)), and therefore somewhat diffi- 
cult CO estimate accurately, but non-resonant quisi-linear theory would sug- 
gest 


3V 


2 




(45) 


where Is the wave potential. This has been estimated by Dum and Duprees 


2 1/2 

<|e > < 0.2 


(46) 


With kj^ R S5/V, kj/-’j^ <0.1 [Kindel and Kennel, 1971] and T^j « 1 heV, (46) 
indicates an RMS wave amplitude of tens of mV/m (as soon by auroral satel- 
lites; Temerin et al., [1981]) and (45) yields a slowing-down time of (1- 
10)0"^. This is a lower limit on the e-folding time for ion heating in Tj^, 
which is driven by electron energy losses. It is a very generous lower limit. 
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not netually am'ruaelUHl in iaet; it eorroopondo to oloiiar naivo eotlmateo of 
Y, « U.l n or of Ion UeatlnB rateo whtoh do not oonoidop the Uraltationn on 
energy tranofor {itoeeoneH ariaing from glt>bal conotdei‘attt>n«t The only point 
In making uheoe eBtimateo io to ohov? that the ratefj oi ttainsfer trora one 
partielc opecloa to another, or to waveo, are far greater than needed, a point 
also eraphaslaed by I’alniadeooo et al. 11974 j. 

Shortly we will make estimates of the actual transfer rates on auroral 
field lines, getting much smaller values. In this conneetlon it is important 
to recall that, with ion resonance broadening as a dominant saturation roeeha“ 

a 

nlom, the ion diffusion Goeffieient I io not ottongly coupled to 

very small values of U are consistent with the Dum-bupree estimates of as 

in (40). 


Palmadcsso et al. 1 1974] hove pointed out that ion heating reduces the 
growth rate Yj^» thus taking the system toward marginal stability, that Is, 
Y “ bat with reduced values of both y and bet us make a quanti- 

tative estimate of this effect. Consider a plasma which is imstoble to wave 
growth, but in which there are no waves to begin with. The condition of 
Instability requires > Vq, vdioro the critical drift velocity Vq is simply 
the minimum drift velocity (over k^ , kj^) . KJ.ndel and Kennel 11971] find the 
minimum at kj^ V » n, kj/kj^ « 0.1 for an isotropic plasma with T - T^. With 
geometric resonance broadening added, one finds for this isotropic case 
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(47) 


The geometric resonance broadening correction is normally 20/o - 30^ of the 
Kindel-Kennel term. The frequency u depends on (see (39), (40); inserting 
this dependence in (47) and differentiating yields the estimate 
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(«) 




(oorao unimportant terrao have been omitted) • 

Wc can find the time rate of change of V^, by Unowlng the ion heating 
rate, which Is roughly MU per particle as one finds from the Boltcmann equa- 


tion (10) by multiplying by 1/2 M and Integrating # 

2 2 ^2 

D e •yji /k^ « yVj^/2, the ion heating rote is 


With 


1^/1^ • If 


< 49 ) 


and from (48) 


V,/Vo" Y I^/T„ « Y 


(50) 


Using (40) - (42), we estimate (setting all temperatures equal) 


1/2 ^ 

Y « Tf ' (Vq - Vg) 


(31) 


So equations (50) - (51) lead to a differential equation for V^; 


, 1/2 ^ ^1 




^C> 


(52) 


Now consider the scenarios where Vjj is constant, which is appropriate to 
many laboratory experiments and possibly some auroral eases. Equation (52) is 
readily solved, showing an exponential approach to equilibrium (V^ » 0). The 
characteristic time is 
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( 53 ) 


A 1 / 1? ^ 1 ^ U I "1 

T « - ^ 1 « 20 n ^ 

V 

e 

wliiirc the value toe t otorao from the .'Otiraate « 0.2 f,lvon by Kindel and 

Kennel (1971) for T„ « T. It followo that the plasma can deviate cubotantial” 
ly from marginal stability only for times ^ T " 0.2 - 2 oec> following oomo 
initial, unstable disturbance. 


Actually, it is quite unrealiotlc to suppose that Vp remains constant as 
Vq changes. The same effect - ion heating - which tends to increase also 
tends to reduce the drift velocity, as we now estimate. Begin with an over- 
simplified picture of the change in the quaol-noutrality condition, in which 
the heating of electrons in is ignored. A change in ion temperature 
must be accompanied by a change in the electrostatic potential <}), in order 


that quasi— neutrality can be maintained! 


\ (Tx '^V ^ ^ \ 'I* + 

Originally this equation was satisfied with 6Tj^ and 6<}> equal to zero. If the 
distribution functions are Maxwellian, and expressed in terms of constants of 
the motion as in (12) - vU), the first-order expansion of (54) is, very 
roughly, 
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(55) 


All we hope to do here is to get the signs right j each term in (55) might be 
multiplied by a numerical coefficient somewhat different from unity. With 
a Ng, Tj^ » Tg, (55) yields 
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(36) 
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The change in the electron drift velocity is directly obtained, because 1/2 ra 
2 

Vj. + pB - e<}) is constant and the magnetic moment p is unchanged; 


(SV,. 


« + 


c <S<{> 
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(57) 


This shows that decreases with an increase in Tj^, but it is not straight- 
forward to incorporate (57) into the differential equation (52) for Vq. The 
reason is that information concerning the changes 6^, 6Tj^ is not transmitted 
instantaneously to the whole auroral field line; instead it must propagate 
along the line at a finite velocity. Clearly the characteristic time for 
changing a parameter like which characterizes the global ion distribution 
function is the inertial time A/V, which is several seconds and much larger 

A 

than the time it takes to reach marginal equilibrium,!. Thus in (52) 
should not be evaluated at the same time t as is V^, but rather at a retarded 
time. In spite of this iechnical complication, we can appreciate from (52) 
that it is possible for to go from positive to negative in a finite time if 
Vq is decreasing, rather than simply going asymptotically to zero as in the 
constant Vq case. If this does happen, the growth rate becomes negative, and 
turbulence originally present is damped away. At this point, ion heating 
Stops, to be replaced by ion cooling at a rate Tj^/Tj^ “ - A/V as new ions are 
injected. After roughly an inertial time, the ions will be cool enough 
in Tj^ to be unstable once again, and more turbulence is generated. We have, 
then, a relaxation oscillator to generate auroral pulsations, with a charac- 
teristic time constant V/A 6 - 10 sec. 
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A complete discussion of this mechanism - including damping - will be 
given elsewhere# Here we only comment that 1) whether oscillations actually 
take place or not depends sensitively on parameters and initial conditions j 
and auroras are not required to pulsate} ?.) our mechanism and its time scale 
are quite different from a recent proposal for auroral flickering [Silevitch, 
1980], but not necessarily Incompatible with the Coroniti-Kennel [1970] 
pulsation picture. 

We conclude this section with some estimates of plasma parameters when 
marginal stability has set in. It has already been observed in Section I that 
the ion heating in Tj^ during the inertial (« quarter-bounce) time cannot 
exceed the initial ion energy, since that is what drives the heating (by 
making a parallel potential drop and an unstable current). This led to the 
estimate D < V^/(AjO 10^^ cm'^ sec”^, and Cvia (7) or (44)) to Y “ V/A. When 
D reaches this value, ion heating and cooling (by injection of more ions) are 
roughly in balance. 

Then (51) yields 

(Vd - V^)/V^ » v/(jinv^ « 0.02 (58) 

These estimates are crude at best, but they suggest that the drift velocity 
and the critical velocity are closely tied together, and a relatively small 
change in either one can turn off the turbulence. If, for example, Vq = 2 x 
10 cm sec then a change in of 2% of this value, or - 6 Vq » 4 x 10 cm 
sec“'^, reduces wave growth to zero. Such a change might be produced by an 
upward fluctuation in D from its saturation value of “ V^/(4A), producing a 
net ion heating of roughly 
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6Tj_ « 2M / ^ D 


(59) 


1 2 

When (59) is used in (57), one finds an upward fluctuation of »< 4 x cm* 
sec""^ will turn off wave turbulence. 

In the next Section, these crude estimates are made somewhat more quanti- 
tative by computer calculations. 
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IV. DETAILED NUMERICAL ESTIMATES 


To do better than the rough guesses of Section III, we need to know two 
quantities with some accuracy! 1) the relationship between the total poten- 
tial drop <|> and the perpendicular ion temperature Tj^; 2) the relationship 
between a given D(s) and the change in Tj^ for a single pass of ions between 
equator and ionosphere, as given in equations (28) - (29). Because explicit 
diffusion is nearly equivalent to a renormalization of (f) and Tj^ in the ion 
distribution function, the first relation can be found by a straightforward 
modification of the Ghiu-Schulz (1978) calculations; the second is found by 
numerical integration of the differential equations. These calculations are 
reported here, and they show that there is a maximum diffusion coefficient 
above which there is no potential drop. The value of Dj^^ is roughly consis- 
tent with our earlier guesses. 

To achieve the above two objectives, we need to obtain a direct non- 
linear differential equation for 0j^ in the variable B since the Chiu-Schulz 
calculations are best expressed in the magnetic field variable. With a 
straightforward, but slightly different set of manipulations from that follow- 
ing (28) and (29), we obtain 


dB 


(B0j^) = 0 



(60) 


where (ds/dB) for a dipole field can be numerically expressed as a function of 
B to very high accuracy. For a given K^), the heated ion temperature Tj^ = 
(M/2 0 ) can be found by solving (60) with a highly accurate Runge-Kutta 
routine. But <j)(B) must be obtained through the Chiu-Schulz quasi-neutrality 
solution, thus ideally 0j^(B) and <j>(B) should be obtained by simultaneous 
solution of (60) and quasi-neutrality. However, the accuracy to which (60) is 
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derived (first order in D) does not justify such an elaborate and expensive 
procedure. We have decoupled (60) from the quasi-neutrality calculation by 
assuming for <|> in Ujj of (60) 

<|.(B) => A<|. (61) 


where A<j> is the initial potential drop. Since U|| is dominated by the magneto- 
spheric ion thermal energy (~ 10 keV), the use of (61), in which eA^) ~ 1 keV, 
is not a significant source of error. The boundary conditions of (uO) for 
downgoing and upgoing ions must be carefully distinguished. For a downgoing 
ion, the plasmasheet ion temperature determines the boundary value 0 j^q, but 
for an upgoing ion of magnetospheric origin, the boundary temperature at 
s“i is determined in its previous downward trajectory. Theoretically, a 
mirroring ion can keep on heating up to higher and higher temperatures by 
virtue of its bounce raotien; however, we note that the electrostatic potential 
is supported on the bounce time scale [Chiu and Cornwall, 1980] so there is no 
reason to assume magnetospheric ions to be trapped for more than a couple of 
bounces. The diffusion coefficient D in (60) is in general an unknown func- 
tion of s or B. We adopt a simple model such that D is a nonzero constant 
throughout the field line. The constant value of D is varied to determine the 
effect of cyclotron turbulence on the potential (|). Constant turbulence along 


the whole field line is not necessarily realistic, but it will serve for this 

2 

initial inquiry. Note that the 0j^ term in (60) is always positive in the 

direction of particle motion and 0j^ “ nature of the solution, by 

2 

virtue of the negative sign preceding the 0^^ term, is that 0j^ decreases from 


a starting boundary value 


%0 


direction 


motion: 


thus, Tj^ increases from the boundary value. However, the decrease in 0^^ can- 
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not be unbounded, for if 1) is too large, can go through a aero at some B. 
Such a solution is unphysical and serves as an upper bound on D, although this 
upper bound is not yet because the effects of heat diffusion upon needs 
to be determined by solving the quasinoutrality portion of the problem, for 
reasonable values of the diffusion constant (D < 10 cm /sec ), Oj^ usually 
does not go through aero; hence may be determined by whether a quasineu- 
tral solution of particle and field equilibrium can be supported. 

To motivate the detailed calculations of ion heating effects on the 
auroral acceleration potential, let us illustrate the basic dependence of 
magnotospheric ion density distribution as function of potential drop and 
temperature anisotropy. Figure 2 shows the densities n^,j^ of down-going and 
up-going magnetospheric ions at the midpoint of the L°8.4 field line. The 
equi-density contours are shown as functions of the parallel potential drop 
and the perpendicular temperature; since the parallel temperature is held 
fixed at 3 keV, the variation of perpendicular temperature is equivalent to 
variation of anisotropy. We note that, for fixed potential drop, the higher 
the anisotropy (Tj^/T^) the lower the density of magnetospheric ions at the 
midpoint of the field line because the ion mirror shrinks towards Che equator, 
as qualitatively expressed in (60). Similarly, as we would expect, magneto- 
spheric ions are repelled towards the equator by higher auroral parallel 
potential drops which accelerate electrons downwards. Now quasi-neutrality 
requires a certain number of ions to be present to balance against the density 
of electrons; thus, for a given electron density distribution, heating of ions 
would require the potentil drop to decrease accordingly in order to maintain 
the ion density along one of the equi-density contours of Figure 2. This is 
the basic effect investigated in detail below, although the results of Figure 
2 is not directly applicable because the effects of a potential decrease upon 
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the distribution of oloetrons must bo taken into %account also. To illustrate 
this point further, we have Investlgatod the dependence of the potential drop 
upon the ion anisotropy by obtaining solutions of a aUu-Behuli; kinetic model 
for L“8.4 with the following fixed parameters: Tj^_^ ° 7 keV, T^_ “ 3 keV, T|j « 
3 keV, n^j_ ° O.b enr^, ng^, ° 70 cm"^ and njj^, « 5 cm“^# Figure 3 shows the 
result of a search of the acceptable potential drop in the above kinetic model 
as the magnetospheric ion anisotropy is varied by varying the constant 
temperature Tj^ while all ocher parameters are held fixed. We note that the 
potential drop maintainable by the electron temperature anisotropy does 
decrease as the ion anisotropy approaches Che electron value. The stippled 
area shows the uncertainty in our solution search due to our not using a 
special high-accuracy computational routine, which is expensive to run. 


The illustrations above are based on arbitrary variations of the constant 
ion perpendicular temperature; therefore, they cannot be directly related to 
the effects of ion eyelotrou heating. To do so, we need to solve (60) to 
obtain the perpendicular ion temperature in the form of (M/20j^) as functions 
of the heat diffusion coefficient D and magnetic field ratio B/By. The solu- 
tions of (60) for U « 10^^ cm^/sec^ for downgoing (arrows pointing away from 
the equator, B/Bq “ 1) and upgoing (arrows pointing towards the equator) ions 


are shown in Figure 4. The solutions are obtained for an initial magneto- 
spheric ion distribution with Tj^ = 3 keV and I'n =» 1.5 keV, which 
yields (H/20j^) given by the dashed curve (with 1)°0) in Figure 4. The initial 
state for Che upgoing ions is assumed to be the heated state at b^^/By. Since 
Che heating effect is prominent primarily in the upgoing magnetospheric ions, 
only about one half of the ion density distribution is severely affected by 
the poCentlal-drop-reducing effects discussed in Figures 2 and 3. Thus, we 
would expect the potential drop to vary fairly slowly with increasing D until 
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tlie solution of (OU) is sufficiently close to Oj^ <=> 0, near which no quasi- 
neutral solution can be obtained. 

An illustration of the coraplcte scenario described in this paper is 
obtained by combining the solution of (60), such as shown on Figure 4, with a 
simultaneous search for the consistent auroral potential drop as we vary the 
ion diffusion coefficient D. The initial conditions for the computation are 
the same as that of Figure 3 except that the lu Llal potential drop is set at 
2.20 keV (to increase the sensitivity of the quasi-neutrality solution to 
parameter changes). The dependence of the auroral potential drop on the ion 
diffusion coefficient D in the range (10^^ - 3 x lO^'^) cm^/sec^ is shown on 
Figure 5. As we have expected, the potential drop decreases slowly with in- 
creasing D until y ~ 10^^ cm^/sec^, at which point the upgoing ion heating 
effects (Figuie 4) become prominent. Above D = 3 x 10^*^ cm^/sec^, the solu- 
tion to (60) breaks down; while above ~ 2.5 x cm^/sec^ the quasineutral 
solution is no longer smoothly varying, so we leave it indeterminate. Thus, 
wo estimate to be ~ 3 x 10^^ cm^/sec^. Again, as in Figure 3, we have 
not used the high-accuracy mode of the Chlu-Schulz model in order to save 
computation time; this results in some uncertainty indicated by the stippling. 

by investigating the properties of the Chiu-Schulz kinetic model of the 
auroral potential drop, we have demonstrated that ion heating leads to a 
decrease of the auroral potential drop, with concomitant effects as discussed 
elsewhere. 
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V. 


SUmiARY At<D CONCLUSIONS 


We have given a ocmi~quantitative picture of Bclf-eonsistent ion 
cyclotron turbulence on auroral field lines, with the following major 
conclusions : 

1. Unless plasmasheet ions are very much more enorgotie than the elec- 
trons, anomalous resistivity is not a large contributor to parallel electro- 
static potential drops; supporting the kind of potential drop actually 
observed requires too much dissipation of energy to be provided by 1 \put from 
the plasmasheet. 

2. Nonetheless, wave turbulence can be present; the ion cyclotron 
turbulence levels suggested by the ion-resonance-broadening saturation mecha- 
nism of Dum and Dupree are comparable to those observed on auroral field 
lines. 

3. The diffusion coefficient 0 ana net growth rate y are very much 
smaller than estimates based solely on local plasma properties; instead they 

3 

are scaled with global paramf^ters; D ~ V /A, y ~ V/A where V is a typical ion 
velocity and A the length of the portion of the field line on which there is 
turbulence. Crude estimates of these effects are supported by calculations 
based on the Chiu-Schulz quasi-neutrality model, with augmented because of 
diffusive heating of ions in Vj^. 

4. The underlying reason for this limit on U and y is that heating in 
Tj^ must be supported by ion kinetic energy, so the effect is one of transfer 

of parallel energy to perpendicular energy (with losses to electron heating 
and wave generation). This necessarily inhibits formation of a potential drop 
as the ions gain anisotropy and mirror higher along the line; eventually the 
potential is too small to maintain sufficient electron current to drive the 
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inacubility. Thus turbulente decroasoQ the auroral pocoutinl Uropi which is 


opposite to the effect of anomalous resistivity. 

S. Depending on initial and boundary conditions, an unstable auroral 
system may evolve in loss than a few seconds to marginally stable equilibrium, 
with turbulence levels, l), and Y as estimated above. Or it may undergo relax- 
ation oscillations as ion heating increases the critical drift velocity V^. at 
the same time it reduces the actual drift velocity Vy (by lowering the poten- 
tial drop). There are generally two characteristic times for these relaxation 
oscillations, the larger of which is the ion inertial time i/V which is the 
time it takes to restore an unstable situation from one in which Vy - Vq has 
turned negative. This time is of order 3-lU seconds, which is observed to be 
the dominant period of auroral pulsations. In another v^ork we will consider 
the pulsation problem in more detail, contrasting our mechanism with earlier 
proposals [Coroniti and Kennel, I970j dilevitch, 1980]. 
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Appendix 


In this appendix we give an alteniacivo development of an approximate 
Qolutlon to (lU) which yields eoacntlally tlie same results as in Section II. 
Since this dovolopraont is based upon different arguments, we offer it hero for 
alternative insight. 

Substitution of the chosen form (2A) into ( 10) yields 
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where *f(s) s ri(s) + (2e/M) 0 |j ipCs) and 6^ is tacitly assumed constant. How 
wo take the n*^*' moment of (Al) in Vj^. Since / dVj^ and / dVj^ Vj_*'“’^^f 

cannot maintain a unique ratio to each other for all n, we conclude that (dA) 
satisfies (lU) only if 
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r au 9 

CA3) 

Elimination of the D-assoeiated terms in (A2) and (A3) yields 

30. . 

^ “■3r--^Blf <®l -^“**^171 f “0 (A4) 

{ Vj ||+ 4D0j f « 0 (A5) 
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In (A4), the spatial and velocity operators on f have been separated into 
product form; therefore, the curly bracket i.iust itself be zero. The solution 
to (A4) is just (31). liquation (A3) is the equivalent of (29). Its solution 
is non-trlvlal because the operator mixes space and velocity space func- 
tions. The solution ansatz to (29) is therefore equivalent to replacing the 
velocity space variable Vj by the spatial variable U|j , in which case the curly 
bracket in (A5) is itself zero. 


40 



Flfiuce Captions 


Figure 1; 
Figure 2: 

Figure 3: 

Figure 4; 

Figure 5. 


The coordinate system and dipole geometry ossumod in this paper. 
Constant density contours of upgoing and downgoing magnetospheric 
ions at the midpoint of the L°8.4 field line. 

Auroral potential drop as function of ion temperature anisotropy in 
the Chiu“Schulz kinetic model. 

Perpendicular ion temperatures on the L“8.4 field line as 

Indicated by the ratio of the local magnetic field B to the equa- 
torial magnetic field Bq. 

Auroral potential drop as function of ion diffusion coefficient. 
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